function PM= Fangbo1(N,rc)%N表示方波的周期数量，rc控制上升沿下降沿的占比，rc越大，占空比越小
t1=0:0.001:1; %周期为1的点，作为时间序列
vt=1-exp(-(t1)/rc);
vt1=exp(-t1/rc);
Y = circshift(vt1,500);
f0=1;         %调制信号频率
pm_t=0.5.*square(2.*pi.*f0*t1)+0.5;%生成方波
Z=pm_t+(Y+vt)-1;
for i=0:N
    y=mod(i,1000);
    if y==1
        Z(i)=0;
    end
end
t=0:0.001:N;%时间序列
pm=0.5.*square(2.*pi.*f0*t)+0.5;  %生成方波
for i=1:N
    for j=1:1000
        pm(i*1000+j)=Z(j);
    end
end
PM=pm;
end